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Abstract 

Tunneiing conductance spectra between a normai metai / c?-wave supercon- 
ductor junction under the presence of buik impurities in the superconductor 
are studied. The quasiclassical theory has been apphed to calculate the spa- 
tial variation of the pair potential and the effect of impurity scattering has 
been introduced by i-matrix approximation. The magnitude of a subdomi- 
nant s-wave component at the interface is shown to robust against the impu- 
rity scattering while that for a subdominant d xy -w&ve component is largely 
suppressed with the increase of the impurity scattering rate. The zero-bias 
conductance peak due to the zero-energy Andreev bound states is significantly 
broadened for the case of Born limit impurity compared with that of unitary 
limit impurity. 



Recent extensive experimental and theoretical studies have revealed that the pair poten- 
tials of high-T c superconductors have d x 2_ y 2-wave symmetryErill. One of the most remark- 
able differences of 0^2 _ 3/ 2- wave superconductors from conventional s-wave superconductors 
is the presence of internal phase of the pair potential. As for tunneling phenomena, the ap- 
pearance of zero-bias conductance peak (ZBCP) in tunneling conductance at [110] surface of 
d x 2_ y 2-wave superconductors reflects the sign change of effective pair potential through the 
reflection of quasiparticle at the surface.^ Qualitative features on the orientational depen- 
dence of the tunneling spectra in high-T c superconductors have been experimentally checked 
by several groups and good consistencies have been obtainediHl!. However, detailed line 
shape of the conductance spectra are not fully understood yet. For example, the ZBCP is 
expected to become ^-functional form for the high-barrier limit case based on the theoretical 
formula, while experimental spectra are largely broadened in usual cases and the origin of 
the broadening is still an open problem. 

Possible explanations for the broadening are the interface roughness and the impurity 
scattering which may inevitably exist in the actual samples. Several existing theories clarified 
the influence of the diffuse scattering at the surface on the surface density of state0~0 and 
on the Josephson current0@. On the other hand, Poenicke, Barash, Bruder and Istyukov 
(PBBI theory) studied the influence of the bulk impurity on the surface density of states 
and Josephson current^. They discovered that the broadening of Andreev bound states due 
to unitary scattering is substantially weaker than that in Born limit scattering. Although 
PBBI theory clarified important aspects of the bulk impurity effect on the charge transport 
effect on the <i-wave superconductor junctions, there still remain several open problems. One 
is the detailed dependence of the conductance spectra on transparency of the junction. In 
PBBI theory, only the junctions with low transparency limit is treated. However, trans- 
parency of the junctions in actual experiments is not restricted to this limit. Actually, it 
has been revealed by the previous theories that the line shape of the tunneling conductance 
strongly depends on the transparency of the junctionifl'i. The other is the effect of the bulk 
impurity on the local inducement of the subdominant pair potentials which break the time 



reversal symmetry. This is important because recent theoretical and experimental study 
have shown the possibility that a subdominant s-wave component or d xy -weive component is 
locally induced near the interface of d x 2_ y 2-wave superconductorS&B. On the other hand, 
the appearance of the predicted splitting of the ZBCP's has been detected only in several 
restricted experiments and conflicting results have been reported thus far@JiH. Therefore 
it is an interesting topic to study the stability of the subdominant pair potential under the 
existence of bulk impurity. 

In this paper, we will calculate the tunneling conductance in normal metal / insulator / 
d x 2_ y 2-wave (n/I/d) superconductor junctions using the quasiclassical formalism of uncon- 
ventional superconductors. The self-energy is calculated in the framework of self-consistent 
i-matrix approximation. We concentrate on two limits, i.e. the limits of the Born scattering 
and the unitary scattering. For [100] oriented junction, the obtained tunneling conductance 
at zero bias voltage is enhanced (suppressed) for high (low) transparency of the junction. 
The degree of the enhancement (suppression) is much more significant for the unitary scat- 
tering case. While for [110] oriented junction, the broadening of the zero bias conductance 
peak (ZBCP) is much more significant for the Born scattering casefi We further study 
the situation where a subdominant s-wave component which breaks the time reversal sym- 
metry is induced near the interface of d-w&ve superconductor. The subdominant s-wave 
component seems to be robust against the presence of the impurity scattering. While the 
subdominant d xy -wave component is reduced drastically with the increase of the magnitude 
of the scattering rate of impurity. 

The model examined here is a two-dimensional n/I/d junction within the quasiclassical 
formalism!!!! where the pair potential has a spatial depends 



lence 
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Here 9 is the angle of quasiparticle trajectory measured from the x axis. If we apply this 
formula to d-wave superconductors including a subdominant s-wave and d xy -wave component 
near the interface, Ar(x, 6) is decomposed into 



A R (x, 9) = A d (x) cos[2{6 - a)} + A s (x) + A d ,(x) sin[2(# - a)] (2) 

where a denotes the angle between the normal to the interface and the x axis of the crystal. 
The insulator located between the normal metal and the superconductor is modeled by a 
5 function. The magnitude of the 5-function denoted as H determines the transparency of 
the junction a^, with cr/v = 4 cos 2 0/{Z 2 + 4 cos 2 9] and Z = 2mHj7\ 1 kp. The effective mass 
m and Fermi momentum kp are assumed to be constant throughout the junction. We solve 
the Eilenberger equation under the existence of bulk impurity in the superconductorEHS. 
The tunneling conductance is obtained thorough the coefficient of the Andreev and normal 
reflect ion&il. In the following, we will concentrate on the normalized tunneling conduc- 
tance a T {eV^ 
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Note that ax{eV) is expressed only by t]r,±{x, 9) just at the boundary (x = 0) where 

T]R t ±(x,9) obeys the following equations^, 

d , 1 



dx ' ihvF cos 9 
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dx ' ihvp cos 9 



-A^E, x, 9+)v% + (x, 9) - A 2 (E, x, 9 + ) + 2A 3 (£, x, 9) Vr , + (x, 9) 

(5) 

-A 2 (E, x, 9„) v 2 R _(x, 9) - A^E, x, 6 J) + 2A 3 (£, x, 9) m ^{x, 9) 



(6) 

A±(E, x, 9) = A R (x, 9 + ) — a±(E, x) + ia 2 (E, x), 
A 2 (E, x, 9) = A* R (x, 9+) + ai{E, x) + ia 2 (E, x), 

A 3 (E,x)=E-a 3 (E,x) (7) 
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with E = eV, vf = kp/m, 9 + = 9 and 6- = 71 — 6, and A R (x, 9+) [A R {x, 9 J)] is the effective 
pair potential felt by an electron [a hole] like quasiparticle with an electron injection from 
the left normal metal. The quasiparticle energy E is measured from the Fermi energy. 
The spatial dependence of the pair potentials are determined by the following equation 



AJx) 



9sk B Tj2^- f ^ d9'{[g R (9',x)] 12 - [g + R {9' , x)] l2 } 

tT 2lT J-n/2 



(8) 
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A d ,(x) = 9d ,k B Tj2^- I*'* d9'sm[2(9' - a)]{[g R (9' , x)} 12 - , x)]^} (10) 
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with dimensionless inter-electron coupling of the s-wave g s , d^-wave ga> and d x 2_ y 2-w&ve 
gd, respectively. The self-energy a>i(E,x) {i = 1,3) describes the impurity scattering and is 
given by 



a,i(iuj n ,x) 



h <gj(9,x)> 
2t 1-cr 



1- i^E* <<?;(M) >< 



(12) 



9r(9, x) = gi(9, x)?x + g 2 (9, x)f 2 + g 3 (9, x)f 3 

with uo n = 2irkBT(n + 1/2). In the above, fj denotes the Pauli matrix and < . . . > means 
the average over the Fermi surface. Here h/(2r) denotes the normal scattering rate, while 
a measures the strength of a simple impurity potential. The spatial dependence of the 
quasiclassical Green's function g R (9, x) is determined by0 



9r{9, x) = U R (9, x, 0)g R (6, 0)U^(9, x, 0) 



(13) 
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U R (9,x,0), 



(14) 



with Ur(9, 0, 0) = 1. In the actual numerical calculations, r]R t ±(x, 6) is calculated from Eqs. 
(|) to (|6|). Since §r(9,0) is expressed by r]R,±(9,0), gji(9,x) is obtained using Eqs. ( |T3| ) to 
(0). Subsequently, the spatial dependence of the pair potentials Ad(x), A^(x) and A s (x) 
are calculated by Eqs. (||) to (|ITD . The self-energies ai(iu n ,x) are calculated by Eq. (^2]). 
To get self-consistently determined pair potentials and self-energies due to the impurity 
scattering this process is repeated until enough convergence is obtained. In the following 
calculation, the two important parameters characterizing the impurity scattering are 2t ^ - 
and a where A c0 denotes the magnitude of the bulk pair potential without the existence of 
impurity scattering. The magnitude of A depends on 2r ^ - and a. 

First, let us consider the case where neither A s (x) nor A^(x) present. One of the typical 
example is the [100] oriented junction, i.e., a = case, where both the transmitted electron 
like and hole like quasiparticle feel the same sign of the pair potentials. For high transparent 
junctions i.e. Z = [see Fig. 1(a)], only the Andreev scattering takes place at the interface 
and the resulting conductance shows zero-bias enhancement due to the Andreev reflction 
without impurity. With the introduction of the impurity scattering the height of the peak 
is suppressed. The degree of this suppression is much more significant for unitary scattering 
case [see curve c in Fig. 1 (a)]. On the other hand, for larger magnitude of Z [ see Fig. 1(b)], 
ar{eV) has a gap like structure without impurity. The magnitude of ot(0) is enhanced with 
the increase of Z. The degree of this enhancement is much more significant in the unitary 
scattering case. In order to understand these line shapes of the tunneling conductances, we 
will look at eV = 0, where t]r } ±(0,9) = ±i is satisfied without impurity scattering. With 
the introduction of impurity potential, the magnitude of ^r,±(0, 6) is suppressed and we can 
denote rjR t ±(0,9) = ±iexp(— 70). Here, 70 is a function with positive value and enhances 
with the increase of the magnitude of 03(0, 0), i.e., the magnitude of n I . For a = 0, since 

ZTLA c Q 

rj R (0,9 + ) = ±iexp(— 7 ) and 77^(0, 6L) = ±zexp(— 7 ) are satisfied, the resulting <7s(0, 6) 
is [1 + exp(— 270)] and [1 — exp(— 27 )]/[l + exp(— 270)] for Z = and sufficiently larger 
magnitude of Z, respectively. Since the magnitude of a 3 (0, 0) is enhanced for the unitary 
scattering case, the resulting 70 also enhances. Consequently, the magnitude of &s(0,9) is 
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much more enhanced for the Born scattering case for Z = 0, while the opposite situation 
happens for sufficiently larger value of Z. 

Next, let us look at the a = 7r/4 case, where ar{eV) has a zero bias conductance peak 
(ZBCP) with small magnitude of For Z = 0, almost the same result as that in Fig. 

1(a) is obtained. This is because well defined bound states do not exist in the absence of 
the barrier potential. With the increase of the magnitude of Z, the ZBCP shows up due to 
Andreev bound state formation without impurityEHi. As shown in Fig. 2, the height of the 
ZBCP is suppressed with the introduction of the impurity scattering. In the Born limit (see 
curve b in Fig. 2), the reduction of the amplitude of the peak height and the broadening of 
the peak width is most significant. To understand this feature, here, we focus on eV = 0, 
where t] R:+ (0,9) = iexp(— 70) and r] R _(0,9) = — iexp(— 70) are satisfied. The resulting 
<xs(0, 6) is [1 + exp(— 27 )]/[l — exp(— 270)] for sufficiently large magnitude of Z. It should 
be remarked that the magnitude of 03(0, 0) is much more enhanced for Born scattering case 
due to the formation Andreev bound state at the interface^ and the resulting 70 is much 
more enhanced. Consequently, the amplitude of &s(0,9), and or(0) in the unitary limit is 
larger than that in the Born limit. 

In Fig. 3, the width of ZBCP W is plotted as a function of scattering rate fr/(2rA c0 ) 
in the Born limit. For small value of Z, W is insensitive with the increase of %/ {2tA c q). 
However, with the increase of the amplitude of Z, W becomes a monotonically increasing 
function of h/(2rA c0 ). For sufficiently larger value of Z, W is roughly proportional to the 



inverse of y^,/(2rA c0 ) as predicted by PBBI theory^. 

Under the presence of the ZBCP, since the quasiparticle density of states at the zero 
energy near the interface are enhanced, a subdominant s-wave (d xy -w&ve) component of 
the pair potential A s (x) (A^(x) ) can be induced near the interface, when a finite s-wave 
(d^-wave) pairing interaction strength exists, even though the bulk symmetry remains pure 
e^2_ y2 -wave0ii~@. Since the phase difference of the d x 2_ y 2-w&ve and s-wave components 
is not always a multiple of n, the mixed state breaks the time-reversal symmetrySS. It is 
an interesting problem to clarify the stability of this broken time reversal symmetry state 



(BTRSS). Here, we study two sub dominant s-wave component is induced near 

the interface, i.e., A s (x) 7^ 0, A d >(x) = [see Fig. 4(a)], (b) a subdominant d xy -wave 
component is induced near the interface, i.e., A d >(x) 7^ 0, A s (x) = [see Fig. 4(b)]. Here, 
the transition temperature of the s-wave component alone T s and that of d xy -wave component 
T d i are chosen as T s = 0.3T d and T d > = 0.3T d , where T d is the transition temperature of d x 2_ y 2 
pair potential without any other components. The transition temperature T s , T d i and T d 
directly correspond to the magnitude of attractive inter-electron coupling g s , g d > and g d , 
respectively. In Fig. 4, the relative magnitude of the subdominant s-wave or d xy component 
at the interface, i.e., | ImagA s (0) | /A or | ImagA^(O) | /A is plotted as a function of 
h/(2rA c0 ) both in the Born and unitary limit. The magnitude of | ImagA s (0) | /A is 
insensitive with the change of the value of h/(2rA c o), while that for | ImagArf/(0) | /Ao 
is reduced with the increase of the magnitude of 2r ^ - . The resulting ax(eV) is plotted 
in Fig. 5 for 2r ^ - =0.1 both in the Born limit and unitary scattering case. In the case 
where subdominant s-wave component is induced at the interface, the resulting aT(eV) has 
a ZBCP splitting independent of the introduction of the impurity scattering. On the other 
hand, the ZBCP is recovered with the introduction of impurity scattering when subdominant 
pairing is c^-wave (see Fig. 5(b)). In this case, the onset of the ZBCP splitting depends 
on the magnitude of 2r ^ - . From these behaviors, we can distinguish the symmetry of the 
subdominant pair potential. 

In this paper, tunneling conductance between normal metal / <i-wave superconductor 
junctions is studied under the presence of bulk impurities in superconductors. We use 
quasiclassical theory and include impurity scattering by t-matrix approximation, and the 
spatial variation of the pair potentials and self-energy are determined self-consistently. For 
[100] oriented junction, obtained tunneling conductance at zero bias voltage is enhanced 
(suppressed) for high (low) transparency of the junction. The degree of the enhancement 
(suppression) is much more significant for the unitary scattering case. While for [110] ori- 
ented low transparent junctions, the broadening of the zero bias conductance peak (ZBCP) 
is much more significant for the Born scattering case@. We further study the situation 
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where a subdominant component which breaks the time reversal symmetry is induced near 
the interface of d-wave superconductoillllH^il. The subdominant s-wave component is ro- 
bust against the introduction of the impurity scattering. While the subdominant d xy -wave 
component is reduced drastically with the increase of the magnitude of the scattering rate 
of impurity. We will comment on the relation between present results with experimental 
ones. In the cases of high-T c superconductors, impurities intrinsically exist in real samples. 
For example, the enhancement of zero-bias conductance is reported in c-axis observation of 
STM experiments^. We believe that the observed features (especially the enhancement of 
conductance) in this experiment are consistent with the result for the Born limit case. On 
the other hand, in the cases of in-plane tunneling, only a few papers reported the splitting 
of ZBCP in YBCO junctionals and the others report the absence of the splitting (for 
example,!!). If the appearance and disappearance of the splitting is simply determined by 
the Born impurity concentration in samples, the promising subdominant component is con- 
cluded to be d xy -w&ve. However, from the microscopic calculation using the t — J model, the 
subdominant d xy -w&ve is hard to be realized while the subdominant s-wave component is 
much more plausibleBil. In this sence, we should perform the tunneling experiments with 
controlled impurity inclusion in the high-T c cuprates. Also we must reveal the origin of the 
attractive potential which induces subdominant d xy -wave component in order to overcome 
this contradiction. Furthermore, it should be pointed out that the magnetically active inter- 
face also induces the ZBCP splittingEU'El. To resolve the origin of the ZBCP splitting in the 
actual junctions, we must study much more about microscopic electronic properties about 
the interface of high-T c cuprates. 

Throughout this paper, the roughness of the interface is not taken into account. Since 
the roughness of the interface suppresses the magnitude of the zBCpli, the low voltage 
behavior of the o"r(eV) is also expected to be influenced. It is an interesting future problem 
to calculate ar{eV) both under the existence of surface roughness and the bulk impurity. 
In the present paper, the effect of bulk impurity is studied where the spatial dependence of 
t is not taken into account. In the actual sample, it is expected that the value of r near 



the interface is much more smaller than that in the bulk. The influence of Andreev bound 
state on the Josephson effect in <i-wave superconductors is also an interesting issueEHHilEl. 
To resolve the bulk impurity effect on the Josephson current in the junctions with sufficient 
transparency is an interesting future problem since the experimentally accessible junction is 
not always in the tunneling limitH. 

This work was partially supported by a Grant-in-aid for Scientific Research from the 
Ministry of Education, Science, Sports and Culture. 
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Figure Captions 

Fig. I. a T (eV) for n/I/d junction with a — 0. (a) Z — 0, (h)Z = 5. a: clean lmit, 
2r ^ - = 0, b: Born limit, 2r ^ - = 0.1, <r = 0, c: almost unitary limit, 2r ^ - = 0.1, a = 0.99. 
Fig. 2. (T T (e\/) for n/I/d junction with a = n/A and Z = 5. a: clean lmit, 2r ^ - = 0, b: 
Born limit, „ % = 0.1, <j = 0, c: almost unitary limit, | =0.1, <r = 0.99. 
Fig. 3. The width of ZBCP W/A c0 for n/I/d junction with a = n/4 is plotted as a function 
of a: Z = 1, b:Z = 5, and c: Z = 20. 

Fig. 4. The relative amplitude of the pair potential of a subdominant s-wave component 
and d X y-w&ve component for n/I/d junction with a = ix/A and Z = 5 is plotted as a 
function of scattering strength of the impurity 2r ^ - . (a) subdominant s-wave component, 
| ImagA s (0) | /A . (b) subdominant d^-wave component, | ImagA d /(0) | /A . a: Born 
limit, a = 0, b: almost unitary limit, a = 0.99. 

Fig. 5. ar(eV) for n/I/d junction with a = 7r/4 and Z = 5 under the existence of a 
subdominant s-wave component [(a)] and ci^-wave component [(b)]. a: clean lmit, 2r ^ - = 0, 
b: Born limit, 2 4 = 0.1, <r = 0, c: almost unitary limit, 2 4 = 0.1, a = 0.99. 
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Fig. 2 




Fig. 3 



